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Mean-field backward doubly stochastic differential equations (MF-BDSDEs, for short) 
| are introduced and studied. The existence and uniqueness of solutions for MF-BDSDEs is 

established. One probabilistic interpretation for the solutions to a class of nonlocal stochastic 
partial differential equations (SPDEs, for short) is given. A Pontryagin's type maximum 
principle is established for optimal control problem of MF-BDSDEs. Finally, one backward 
linear quadratic problem of mean-field type is discussed to illustrate the direct application 
of above maximum principle. 
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1 Introduction 

Backward doubly stochastic differential equation (BDSDE for short) of the form 

Y t = S+ / f(s,Y s ,Z s )ds + / g(s,Y s ,Z s )dB s - / Z s dW s , < t < T, 
Jt Jt Jt 

was firstly initiated by Pardoux-Peng [22j to give probabilistic interpretation for the solutions of 
a class of semilinear stochastic partial differential equations (SPDEs for short). BDSDEs have 
not only emerged as a natural and convenient tool in the context of SPDEs, see Bally-Matoussi 
[5], Hu-Ren [14], Pardoux-Peng [22], Ren-Lin-Hu [23], Zhang-Zhao ( [26] , [27] ) , but also recently 
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gained interest in other fields as well, especially in relations to the stochastic optimal control 
problems, see Bahlali-Gherbal jlj, Han-Peng-Wu [13], Zhang-Shi [25] , 
McKean-Vlasov stochastic differential equation of the form 

dX{t) =b(X(t),n(t))dt + dW(t), t e [0,T], X{0) = x, (1.1) 

where 



b(X(t),n(t))= / b(X(t,u>),X(t;u/))P(du/)=E[b{t,X(t))] 

Jn ?=-*(*) 

b : R m x R — > H being a (locally) bounded Borel measurable function and fj,(t; •) being the 
probability distribution of the unknown process X(t), was suggested by Kac [16] and firstly 
studied by McKean [20]. So far numerous works has been done on McKean-Vlasov type SDEs 
and applications, see for example, Ahmed |lj, Ahmed-Ding [2], Borkar-Kumar [6], Chan |10| . 
Crisan-Xiong [11] . Kotelenez [171 , Potelenez-Kurtz [18], and so on. It is worthy to point out 
that (jl.ip is a particular case of the following general version, 

X(t) = x + j\( S ,X(s)^<p b ls,X(s),^ =x{s) )ds 

+ J a(s, X(s), E0 CT [s, X(s),^ =x{s) )dW s , (1.2) 

which can be regarded as a natural generalization of classical SDEs. Mathematical mean field 
approaches play a crucial role in diverse areas, such as physics, chemistry, economics, finance 
and games theory, see for example Lasry-Lions |19| . Dawson [T2], Huang-Malhame-Caines [15] . 
In a recent work of Buckdahn-Djehiche-Li-Peng [8], a notion of mean-field backward stochastic 
differential equation (MF-BSDE for short) of the form 

Y t = c + ^ T E7( s , w , w ^y s H,z s H,y s (w , ),^(^))^-^ T z s ^, 

with t G [0, T] was introduced to investigate one special mean-field problem in a purely stochastic 
approach. 

In this paper, we would like to introduce mean-field backward doubly stochastic differential 
equation (MF-BDSDE for short) 

Y t = f(s,Y s ,Z s ,Tf(s,Y s ,Z s ))ds 

rp rp 

+ [ g(s,Y s ,Z s ,T 9 (s,Y s ,Z s ))dB s - [ Z s dV& S) (1.3) 
Jt Jt 

where 

[F 1 (s,Y s ,Z s )}(uj)= [ e\s,u,J,Y s {u),Z s {u),Y s (u'),Z s {J))V{du:'), (1.4) 

with I = f, g. For convenience, we also denote 

Ey (s, Y s , Z S X, Z's)\ ■= r'(s, Y„, Z s ), 
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when there is no abuse of notation. Following the basic ideas in |22| . we firstly discuss the exis- 
tence and uniqueness of solutions for MF-BDSDE (jl,3p . which obviously extends the results in 
both [22J and [9]. It is worthy to point out that MF-BDSDEs not just is a natural generalization 
of BSDEs and MF-BSDEs from the view of mathematics. Our study on them also is motivated 
by the problems in the following two aspects. 

As is well-known to us, the study on stochastic partial differential equations have increasingly 
been a popular issue in recent years. As one kind of them, stochastic partial differential equa- 
tions of McKean-Vlasov type were discussed in [18]. In fact, such equations were obtained as 
continuum limit from empirical distribution of a large number of SDEs, coupled with mean-field 
interaction. We also refer the reader to [11] and [T7] for more details along this. On the other 
hand, we would also like to mention the work of Buckdahn-Li-Peng [9] who studied one kind 
of nonlocal deterministic PDEs. In virtue of the "backward semigroup" method they obtained 
the existence and uniqueness of viscosity solution for nonlocal PDEs via mean-field BSDEs (|1.3|) 
in a Markovian framework and McKean-Vlasov forward equation. Motivated by the above two 
cases, in this paper we will give some discussions on one kind of nonlocal stochastic partial 
differential equations. Since our backward equation here is allowed to dependent on Z°' x °(-), 
therefore the nonlocal SPDEs here is not a direct generalization of deterministic PDEs in [9] to 
the stochastic case. Some additional necessary and essential terms are required in our SPDE 
to meet the general case here, see (|4.3p below. On the other hand, comparing with the case 
in [22] and [9], due to the nonlocal property of SPDEs and the interesting measurability of 
the corresponding solution u(t,x), one fundamental and important term u'(t,x) is required to 
meet such general case, see also Remark 4.2 below. Instead of investigating the limit result for 
such equations, we will study the SPDEs from other aspects. A probabilistic interpretation for 
the solution to such kind of SPDEs is derived by a connection between them and decoupled 
forward-backward doubly differential equations of mean-field type, which extends the results in 
|22] to the mean-field case. 

The second motivation stems from the study of optimal control problem and certain stochas- 
tic differential games problems. Some related works along this have followed two main venues. 
On the one hand, optimal control of mean-field (forward) stochastic differential equations was 
discussed in Andersson-Djehiche [3], Buckdahn-Djehiche-Li [7] and Meyer-Brandis-Oksandal- 
Zhou [21 J where stochastic maximum principle were derived as a necessary condition of the 
optimal control. On the other hand, optimal control problem for backward doubly stochastic 
differential system (or forward-backward system) were spread out in [4], [13], [25] where the 
corresponding linear quadratic problem and nonzero sum stochastic differential games were also 
investigated. Inspired by above two case, it is natural for us to consider the optimal control 
problem for backward doubly stochastic system of mean-field type. Since the terms T l with 
I = f)9 ( see (15. ip ) have more general feature than the corresponding one in [3], [7] and [21], we 
have to introduce E* in our adjoint equation, which is slight different from E'. Note that such 
kind of skill also appears in [23]. On the other hand, since 6* , 9 and I are allowed to depend 
on it(-,u/), some new terms appear and is expressed by E* in our maximum principle (see (15. 9\\ ) 
which is totally different from all the previous literature. We believe that such new features will 
lead to some other interesting things and we hope to explore and study them carefully in the 
future. 

The paper is organized as follows. In Section [2J we will present some preliminary notations 
needed in the whole paper. In Section [31 we consider the existence and uniqueness of solution 
for MF-BDSDE (|1.3j) . In Section HI we give the probabilistic interpretation for the solutions 
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to a class of nonlocal SPDEs by means of MF-BDSDE. In Section we discuss one optimal 
control problem of MF-BDSDE. In Section [61 we investigate one backward doubly stochastic 
LQ problem of mean-field type to show one direct application of result in Section [5j 

2 Preliminaries 

Throughout this paper, let T, P) be a complete probability space on which are defined two 
mutually independent Brownian motion {VFi} f>0 and {Bt} t>0 , with value respectively in 1R d 
and PJ. We denote by 

T t = Jf V^, Vt E [0,T], 
where N is the class of P-null sets of T and 

jrW = a { Wr . <r <t}V N, T^ T = a {B T - B r ; t<r<T}vM. 

In this case, the collection {Tt, t E [0, T]} is neither increasing nor decreasing, while {T^; t E [0, T] } 

is an increasing filtration and ^T^ T ;t E [0,T]| is a decreasing filtration. 

Let {ti 2 , T 2 ,V 2 ) = (0 x fi, T® T, P (g)P) be the completion of the product probability space 
of the above (fl, J 7 , P) with itself, where we define J-f = Tt<S>Tt with i E [0, T] and Tt®T being 
the completion of Tt x Jj. It is worthy of noting that any random variable £ = defined on 
can be extended naturally to fi 2 as £'(w,u/ ) = with (w,u/) E ft 2 . For H = P n , ect, let 
L 1 (Q 2 , J 72 , P 2 ; be the set of random variable £ : ft 2 ^ H which is ^-measurable such that 
E 2 |£l = J Q 2 |£(u/,w)|P(dcj')]P(dw) < oo. For any 7/ E L 1 ^ 2 , J 72 , P 2 ; #), we denote 

E'r/(uv)= / t/(o;,o/)P(cZo/), E*t/(-,o;) = / t/(o/, w)P((ia/). 

Particularly, for example, if 771 (w,u/) = 77(0/) and 772(0;, 0/) = 772(0;), then 

E't/i = / ?7i(o/)P(do/) = E771, E*7/ 2 = / 7? 2 (o;)P(do;) = E772. 

Hence, in what follows, E' and E* will be used when we need to distinguish 0/ from uj, which is 
the case of both u and 0/ appearing at the same time. On the one hand, the well definition of E' 
above gives the precise meaning of and V 9 in (jl.4p . On the other hand, it also indicates that, 
for example, the operator T$ is nonlocal in the sense that the value T$ (s,uj,Y(s,uj), Z(s,uj)) of 
rf (s,Y(s), Z(s)) at uj depends on the whole set 

{{Y{s,uj'),Z{s,uj')) I 0/ E ft}, 

not just on (Y(s, uj), Z(s, uj)). 
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At last, we would like to introduce some spaces of functions required in the sequel. 



S 2 ([0,T];]R n ) = [0,T] xQ^WC 1 , 



(pt is Ft -measurable process 



such that E( sup \tpt\ ) < oo}, 

0<t<T 



M 2 (0,T;R n ) = {ipt : [0, T] x Q R n , 



ipt is JVmeasurable process 



f T 

such that IE / \ipt\ 2 dt < oo}, 
Jo 



L 2 (0,J- T ,P;]R n ) = : [0,T] x ft^]R n , 



£ is J-r-measurable random 



variable such that E|£| 2 < oo}. 

3 The unique solvability of MF-BDSDEs 

In this section, we will discuss the existence and uniqueness of adapted solution for MF-BDSDE 
(jl.3p which is rewritten below (for convenience): 



Y t = £ + jf f(s,Y s ,Z s ,rf(s,Y s ,Z s ))ds 

+ F g(s, Y s , Z s ,T°(s, Y s , Z s ))dB s - F Z s c0 s , (3.1) 
Jt Jt 

with I = f, g. Before it, we make the following assumptions. 

(HI) (i) f G L 2 (0,J- T ,P;IEl n ). / : x [0,T] x R n x R nxd x H fcl H n is measurable 
and for all {t,y,z,j) G [0,T] x R/H-nxd+fc^ ^ f(t,u,y,z,i) is J^-measurable. 

# : ft x [0,T] x R" x ]R nxc! x R k2 -> R mxl is measurable and for all (t,y,z,j) G 
[0,T] x ]R^+«x^+fo ) (t >w ) ^ g(t,uj,y,z,j) is .^-measurable. The map 6 l : [0,T] x 

^2 x R 2n+2nx2d ^ R m j g measurable and for all (t,y,Z,y', z') G [0,7] X ]R 2n+2nXM , 

the map (i, cl>,o/) (-)• G l (t,u),u)' ,y,z,y' ,z') is J-f-measurable on [0, T], I = f,g. 

(ii) / and g satisfy uniformly Lipschitz condition to (y, 2,7), that is, there exist positive 
constants Lj Kj and ay with i = y,z, y' , z', 7, j = 1, 2, 3, 4, such that 

2/1,^1, 71) - f(t,y 2 ,z 2 ,72)\ 

< L y \yx - y 2 \ + -M^i - z 2 | + L 7 |7i - 72I, 

2/1,^1, 71) - 5(*,y2, ^2,72) | 2 

< if 2 |yi - y 2 | 2 + ai^i - z 2 \ 2 + a 2 |7i - 72| 2 , 
\9 f {t,uj,uj',y 1 ,z 1 ,y' 1 ,z' 1 ) - 9 f (t,u,u' ,y 2 , z 2 ,y' 2 , z' 2 )\ 

< L y \yi - y 2 \ + L z \zi - z 2 \ + L y i\y[ - y' 2 \ + L z >\z[ - z' 2 \, 
\9 9 (t,u,uj' 1 y 1 ,z 1 ,y' 1 ,z' 1 ) - 6 9 (t,uj,uj' ,y 2 , z 2 ,y' 2 , z' 2 )\ 2 

< Kl\ Vl - y 2 \ 2 + Kl,\y[ - y' 2 \ 2 + a-z\ Zl - z 2 \ 2 + - z' 2 \ 2 , 
V(t,w,a/) G [0,T] x ft 2 , G H 4n+M ,i = 1,2, 
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and 



T 



•T 



E / \E'6 l (t ,u),u) , )\ ds < oo, E / |Zo(s,w)| ds < oo, l = f,g, 



where 6 l (t,u},u)') = 6 l (t, to, to', 0, 0, 0, 0), /o( s ; w ) = ^o( s > w ) 0, 0, 0). Here we assume 
that a\ + a 2 a 3 + a^a^ < 1. 



Remark 3.1 Under (HI), since we have 



\e\t,ui,ui',Y(t,u),Z(t,u),y,z)\ 
< L(l + \Y(t,u)\ + \Z(t,u)\ + \y\ + \z\), l = f,g, 



\T l (t, Y(t), Z(t))\ < L(l + \Y(t)\ + \Z(t)\ + E|Y(t)| + E|Z(t)|), Z = f,g. 
Likewise, for any (Y 1 (-),Z 1 (-)),(Y 2 (-),Z 2 (-)) G S 2 ([0,T};R n ) x M 2 ([0, T};R nxd ), 



\Tf(t,Y 1 (t),Z 1 (t))-Tf{t,Y 2 {t),Z 2 (t))\ 

< LylY^t) - Y 2 {t)\ + L z \Z x {t) - Z 2 (t)\ 
+L y ,E\Y 1 (t) - Y 2 {t)\ + L z M\Z x (t) - Z 2 (t)\, 
\T^t,Y 1 (t),Z 1 (t)) -TV{t,Y 2 (t),Z 2 (t))\ 2 

< i^|Yi(t) - Y 2 (t)\ 2 + K 2 y ,nYi{t) - Y 2 (t)\ 2 
+a 3 \Z 1 {t) - Z 2 (t)\ 2 + a 4 E\Zx{t) - Z 2 (t)\ 2 . 



The above two estimates will play an interesting role in the next theorem. 

Theorem 3.1 Suppose (HI) holds. Then MF-BDSDE (|3.ip admits a unique adapted solution 



(Y,Z) G S 2 ([0,T];]R n ) x M 2 (0, T; R nxd ). 

Proof. For any (y,z) G S 2 (0,T;U n ) x M 2 (0,T;R nxd ), we consider the following MF- 
BDSDE 



According to Theorem 1.1 in [52], there exists a unique pair of solution (Y, Z) G S* 2 ([0, T]; H n ) x 
M 2 (0,T;lR nxd ). Hence, if we define Q(y,z) = (Y,Z), then 9 maps from S 2 ([0, T]; R n ) x 
M 2 (0,T;lR rixa! ) to itself. We now show that G is contractive. To this end, take any (y l ,z % ) G 
S 2 ([0,T];]R n ) x M 2 (0,T;]R nxa! ), (i = 1,2), and let 



consequently, 




(Y\Z*) = Q{y\z l ). 
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We denote by (Y, Z) = (Y 1 - Y 2 ,Z X - Z 2 ) and (y,z) = (y 1 - y 2 ' ,z x - z 2 ). Using Ito formula to 
e^* | Y t \ 2 we have 

Ee^|F t | 2 + E / e^ s \Z s \ 2 ds + ~E / (5e Ps \Y s \ 2 ds 
Jt Jt 

rT f -T 

= 2E / e^ s Y s f(s)ds + E / e^g^ds, 



where 



f(s) = f(s,Y s \zl,E'[O f (s,Y s \zl,yl',zl')}) 

-f(s, Y 2 , Z 2 , M'[ef(s, Y 2 , Z 2 , yf, zf)}), 

g(s) = g{s,Y},ZlEr[e°(8,Y 8 \zl,y?,z?)]) 

-g(s, Y 2 , Z 2 ,TE>[9°(s, Y 2 , Z 2 , yf , z 2 ')}). 



Hence from (HI), Remark 13. II above and the inequality ab < -a 2 + 5b 2 , we have 

o 

Ee^l 2 + EJ e^ s \Z s \ 2 ds + EJ ^ s \Y s \ 2 ds 

< 21eJ e^[(l+L 7 )|y s |(Ljy s | + L 2 |Z s |) + L 7 |F s |(VE|y s |+L / E|z s |)]ds 

+E J e^ s [K 2 (l + a 2 )|Ys| 2 + (ai + a 2 a 3 )\Z s \ 2 + a 2 K y ,]E\y s \ 2 + a 2 a 4 E|z s | 2 ]ds 

< MiE / e^|Y s | 2 ds + M 2 E / e""|Z,| 2 ds 

+M 3 E / e ps \y s \ 2 ds + M 4 E / eHz„| 2 ds, 
Jt Jt 



where 



Mi = K 2 y {l + a 2 ) + {l + L 1 ){L y + L z -C + -L 2 + L 1 L z , ■ C) 



M 2 = il±^I^i + (a 1 + a 2 a3), 
M 3 = ^L 2 , + a^,; M 4 = + a 2 a 4 . 

After some simple calculations, it is easy to see that 

E f ^l e ^\Y s \ 2 ds + V [ T e^\Z s \ 2 ds 

Jt 1 - M 2 Jt 

< _M±- [ E r^-i^^ + E ( T e^\z s \ 2 ds 

By the assumption imposed on ai, is a contraction on M 2 (0,T;]R n+nxd ), thus there is a 
unique fixed point (Y, Z) G M 2 (0, T; ]R«+™ xd ) which is the solution of (|53]> . Moreover, it is easy 
to check that Y(-) G ,S" 2 (0, T; R n ). The proof is complete. □ 
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Remark 3.2 Note that our result here can fully cover the corresponding results in [9] and |22j . 
In fact, if K y = a\ = a<i = or L 7 = a<i = 0, owr result degenerates respectively to the case in 
[9] and [22] . 

Next we introduce a type of forward doubly stochastic differential equation with mean field type 
as follows, 

Pt = V+ I f(s,P s ,Q s ,r f (s,P s ,Q s ))ds 
Jo 

+ / g(s,P s ,Q s ,T 9 (s,P s ,Q s ))(0 s - [ Q s dJ3 s , (3.2) 
Jo Jo 

where rj is J-o-rneasurable. Note that such kind of equations appear as adjoint equation in the 
optimal control problem below. Here we will transform (|3.2p into the similar form of (|3.ip . If 
we define 

B t = B T - B T _ t , W t = W T - W T -u K = T$ V Tf , 
P t = P T -t, Qt = Qt-u te[0,T], 

then Ft = Ft-u V 1S ^r-measurable, Pt, Qt are J^-measurable, and 

Pt = V+ f{s,P s ,Q s ,rf{s,P s ,Q s ))ds 

+ / 5 (s,P s ,Q s ,rf(s,P s ,Q s ))d^ s - / Q s dB s . (3.3) 

Note that (|3.3p have the similar form as f)3. 1 1) . then by Theorem 3.1 there exists a unique pair 
of (P,Q) solving (|3.3p . and we have 

Theorem 3.2 Suppose (HI) holds. Then (|3.2|) admits a unique adapted solution (P,Q) € 
S 2 ([0,T};H n ) x M 2 (0,r;]R n><d ). 



4 Probabilistic interpretation for a class of nonlocal SPDEs 

The connection between BDSDEs and systems of second-order quasilinear SPDEs was firstly 
observed by Pardoux and Peng [22J, where the probabilistic interpretation for second-order 
SPDEs of parabolic types was derived. Thereafter, this subject has attracted a lot of research, 
such as [5], [H], [23], [26], [21]. This section can be regarded as a continuation of such a theme. 
In other words, we will exploit the above theory of MF-BDSDE in order to provide a probabilistic 
formula for the solution of a class of nonlocal SPDE. 

Given arbitrary xq £ H m , (t, x) € [0, T] x R m being the initial condition, let us consider the 
following forward SDE in R m , 

XY = x + f T b {r, Xl' x )dr + f T a {r, X^dV^r, s > t, (4.1) 
Jt Jt 
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and the backward equation 

Yg ,x = Et[h(X% B ,(X%r Bo y)]+ [ T T{(r,X t r > x ,Y?> x ,Zt.> x )dr- f Z^ x (0 r 

J s J s 



+ f T rf (r, x l fX ,x , Z^dSr, s > t, 

J s 



(4.2) 



where 



T k {r,X l r > x ) = / k{r,Xl> x (u),X^ Xa (uj'))V(duj') =: E'[fc(r, X*'*, (X°' x °)')], 
Jn 



\r, X£ x (w), 17' x (u;), ^ x (w), X^o/), ^(u/), Z?' x °(u'))F(du}') 
=: E'[0 J (r, X*'*, Y^, (X r ^ )' ', {Y r °' x °)' , {Z?> x °)% 
with k = b,a, I = f,g, and 



b : [0, r]xR m xR m 4 R m , a : [0, T] x R m x R m 
^ : [0, T] x R m x R n x R nxd x R m x R n x R nxd 



mxd 



9 : [0,71 xR m xE"xE 



nxd 



R m x R" 



R 



nxd 



> R 

R n , 
R nx ' 



/i : R m x R" 



R n . 



It is known that SDE (|4.1|) has a unique solution if coefficients satisfy linear growth and the 
Lipschitz condition, see [9]. Similarly under suitable assumptions of h, O* and 9 9 , we can also 
obtain naturally the wellposedness of (|4.2|) by means of Theorem 3.1 above. 

Suppose 0f(t, x, x', •, ■, •) and 9 9 (t, x, x' , •, -, ■) satisfy the conditions in Theorem 13.11 uniformly 
for t, x and x', and E|/i(X^ a: , (X^ x °)')\ 2 < oo, it follows from Theorem 13.11 that (|4.2|) admits a 
unique solution (Y t , Z t ) G S 2 ([0, T]; R n ) x M 2 (0, T; R nxd ). 

We now relate MF-BDSDE P~2|) to the following nonlocal SPDE: 



du(t, x) 



(Cu(t, x) + E'6 f (t, x, u(t, x), v«(t, a;) • E'[er(t, x, (X t ' x ° )')] , 



(xr°)>'(M^r)')jv«'(^) -Ev^^^xr )')] v o„ P ) * 

a=X t 

-E'6 9 (t,x,u(t, x), \/u(t, x) ■ E'[cr(t, x, (X^ ,X °Y)], (Xt' X0 )', 



u'(t, (X?> x °)'), [ V n'(t,a) • EV(t,a, (X t ^)')] 



0,x 



it(T, x) 



E>[h(x,(X% x °)')}, 



(4.3) 



where it : R x R m x O — > R n , £u = (Lui, • • • Lu n ) T ', with 



Lm(t, x) 



*r(E / [a(t J x J (X°' a,0 ) , )]E , [<7(t,x,(X t u ' !eo )')] r £> 2 «(t ) x)) 



0,3;o\/\iT n 2„ 



+ V«(i,x)-E>(i,x,(X t °> x °)')], 
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and 

u'(t,(X?> Xo y) = u(t,u;',X?> Xo (u,% t£[0,T]. (4.4) 

We can assert that 

Theorem 4.1 Suppose that b, a, f and g satisfy suitable linear growth and Lipschitz condi- 
tion, h xx (x, (X^ x °)) exists and EE' \h(Xj, x , (X^ x °))\ 2 < oo. Suppose SPDE (|4.3p has a solution 
u(t,x) E C 1,2 (fi x [0,T] x H m ;]R n ). Then, for any given (t,x), u(t,x) has the following inter- 
pretation 



u(t,x) = Y t l > x , (4.5) 



where Y t ' is determined by (|4.1|) and (|4.2|) . Moreover, the solution u(t,x) of (|4.3p is unique 
too. 



Proof Applying Ito's formula to u(t,X t ), we obtain 

u(T,X^ x )-u(t,x) = j\^-(r,X^)+Cu(r,X^)]dr 



+ 



V«(r,*?' x ) • E'[a(r,X^, (X r °^)')]o# r . 



Because u(t,x) satisfies SPDE (|4.3p . it holds that 



w(T, ) — u(t, a;) 



- J tie' (s, Xt> x , u(s, X?), V u(s, Xt> x ) ■ Ef[a( 8 , X?, (X° s > x ° )')], 
(X°^)',u'( S , (X°'*°)'), [v«W) • EV( S ,a, (X°>*°)')] )<fa 

a— A a 

- ^ E'^( S , Xl' x ), V «(«, X*'*) • E'[a( S , Xf'*, (X°>*°)')], 
(X°>*°)>'( S , (X s °'*°)'), [V«'(«,a) ' EV( S) a, (X°>*°)')] 

a=Jf s u 



By the uniqueness of solution for (|4.ip and (|4.2p . it is easy to check that (u(s, xt' x ), \/u(s, xt' x ) 
E'[a(s, X*'*, (X^ )')])) with s E [0, T] is a solution of (JOD- Hence it follows that 



u(t,x) = 1?'*, i e [0,T],x E R m . 

Under the above condition, the solution Y t,x {-) is unique, then the solution u(t, x) of SPDE (|4.3p 
is also unique. □ 

Remark 4.1 (|4.5p can 6e regarded a stochastic Feynman-Kac formula for SPDE (|4.3p . which 
is a useful tool in the study of the property for SPDE. For example, from Theorem 4-1 we know 
that the solution of SPDE (|4.3I) must be unique if it exists. 
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Remark 4.2 Note that the introduction of function v! in (|4.4p coincides with the general setting 
in our discussion. Actually, on the one hand, comparing with the SPDEs in |22| . the term 
u'(t, (X^ ,x °)) is necessary since our SPDE here is nonlocal. On the other hand, comparing 
with the case in [9], here we replace u(t, (X^' xo )') frequently used there with a slight new term 
u'(t, (X^ ,x °y) because of the special measurability of u(t,x). As we know, since Y t,x (s,ui) = 
u(s,w,X*' a: (s,a;)), thus Y°> x °(s,oj') = u(s,ui' , X°' x ° (s,u/)) = u'(s,(Xs' x °)'). Obviously, if 09 
equals to zero, i.e, (|4.3p becomes a deterministic PDE, hence u'(t, (X t ' x °)') will degenerates into 
u(t, (X^ ,x °)') which is the case in [9]. 

Remark 4.3 In our framework, (j4.2|) is allowed to depend on Z 0,x °(-), hence some more neces- 
sary terms representing the nonlocal property of (|4.3|) are needed. This is totally different from 
the discussion in [9]. 



5 An optimal control problem for MF-BDSDEs 

In this section, we would like to consider one optimal control problem for MF-BDSDEs. As 
a necessary condition, we will derive one maximum principle. For the reason of simplicity, we 
assume m = n = d = l = ki = k2 = l- 

Given a convex subset U C for any admissible control v 6 U a d, where 



U, 



ad 



jt> : [0,T] x ->■ U\v is T t -measurable, E J \v t \ 2 dt < +oo| 



and £ e L 2 (0,^ T ,P;E), we consider the following MF-BDSDE: 



Y t v = £+ C \ f {s,Y:,Z v s ,v s )ds+ C T9(s,Y s v ,Z*,v s )dt s - F Z, 
Jt Jt Jt 



(5.1) 



where i = f, g, 

r(s,Y s v ,Z v s ,v s ) 

= [ o i ( Sl ^',y;( w ),2:( u ) 1 » ! H,y:M 1 2:M,» s M)p( ( L'), 

Jn 

and 

9 f :fi 2 x[0,T]xRxRx[/xRxRx[/4R, 
« J :Sl 2 x[0,T]xRxRx[/xRxRx[/4R. 

The control problem is to find an admissible control to minimize over U a( { the cost function of 

JK-)) =E [ T T l (s,Y s v ,Z v s ,v s )ds + m'h(Y^(u;),Y v (u J '))}, (5.2) 
J o 

where 

i(s,wy,i7( w ) 1 ^( u ) ) «,( u ) ) i7(w0,^(^),t',(u/))p((iy) 



ii 



and h : tt 2 x R x R R, 

Z : Sl 2 x[0,r]xRxRx[/xRxRx[/-jR. 

The following is our main assumptions on the above mappings in this section: 

(H2) (i) 9* , 9 9 , I and h are continuous and continuously differentiable with respect to y, y' , z, 
z', v, v', and the derivatives of I and h are allowed to be linear growth. 

(ii) 9* and 9 9 satisfy uniformly Lipschitz condition to (y, z, y' , z', v , v'), that is, there exist 
positive constants Lj Ki and atj with i = y,z,y',z',v,v', j = 3,4, such that 

\9 f (t,uj,uj',y 1 ,zi,y' 1 ,z' 1 ,v 1 ,v' 1 ) - 9 f {t, to, to', y 2 , z 2 , y' 2 , z' 2 ,v 2 , v' 2 )\ 

< L y \yi - y 2 \ + L z \zi - z 2 \ + L y >\y[ - y' 2 \ 
+L z '\z'x - z' 2 \ + L v \v\ — v 2 \ + L v i\v[ - v 2 \, 

\9 9 (t, u], J, y 1 ,z 1 ,y[, z'^vx.v^) - 9 9 {t, u, J , y 2 , z 2 ,y' 2 , z' 2 ,v 2 ,v' 2 )\ 2 

< K 2 y \ yi -y 2 \ 2 + K 2 \y[ - y' 2 \ 2 + K 2 v \v x - v 2 \ 2 

i 7^2 I / l|2 , I i2 , I / / |2 

+iv !; ,|w 1 — v 2 | + a 3 \zi - z 2 \ +a 4: \z 1 -z 2 \, 

V(t,w,a/) G [0,T] x ft 2 , 4 G R 6 ,i = 1,2, 

and 

/•T 

E / |E'6^(t,w,u/)| 2 (i,s < oc, l = f,9, 
Jo 

where 9 l (t, to, to') = 9 l Q (t, to, to', 0, 0, 0, 0, 0, 0), and 03 + 04 < 1. 

Under the above hypotheses, for every v(-) G U a d, by Theorem 13. 11 (|5. 1 1) admits a unique strong 
solution (Y V ,Z V ) G S" 2 (0,T;R) x M 2 (0,T;R), andthe cost functional J is well-defined. 

Suppose that u(-) is an optimal control and (Y(-), Z(-)) is the corresponding optimal tra- 
jectory. Let v (•) be such that u(-) + u(-) G U a d- Since is convex, then for any < e < 1, 
u £ {-) = u(-) + ev(-) is also in U a d- Prom Theorem 13. 1\ we know that state equation (|5.ip has 
a unique solution, denoted by (Y £ (■) , Z £ (■)) corresponding to u £ . Before the main result we 
require to prove some basic results. 

Lemma 5.1 Under assumption (H2), for any t G [0, T], we have 

r T 

E\Y £ - Y t \ 2 < Ce 2 , E J \Z% - Z s \ 2 ds < Ce 2 . 
Proof. Notice that Y t £ - % satisfies the following MF-BDSDE: 

Y £ -Y t = J [Tf(s,Y £ ,Z £ s ,u £ s )-Tf(s,Y s ,Z s ,u s )]ds 

+ ^ V(*, Y £ , Zl, u%) - T°(s, %, Z s ,u s )]dJ3 s 

- J\z £ s -z s )dvfr s . 
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Applying Ito's formula to \Y t £ — Y t \ 2 , we have 

M^\Yf-Y t \ 2 + ^ T \ZI-Z s \ 2 ds^ 
= 2E ^ (Y a e -Y s ,rf(s, Yg ,Z £ s ,u £ s )-rf(s,Y s ,Zs,u s )) ds 

+e J T |r»(a, yf, z a e , <) - r»(«, y , z s , u s )| 2 d s . 

From (H2), we have 

E|y/ - y| 2 + Ey |zf - z s \ 2 ds 

< feiE y |y/ -Y s \ 2 ds + k 2 e 2l E £ \v s \ 2 ds, 



where fcj (i = 1,2) are two constants depending on the parameters in (H2). By Gronwall's 
inequality the Burkholder-Davis-Gundy inequality, we obtain the results. □ 
In the following we make the convention that 

= a(•,a;, W ',y £ (•,a;),Z £ (.,t i ;),n £ (•, W ),y £ (•,a;'),^(•,^)^ £ (•^ , )), 
where oj, u/ € f2, a = 9* , 9 9 , 1. We introduce the variational equation as follows, 

tt = i>t+ ! F 1 (s,£ s ,r] s )ds+ [ Gi(s,^,r Js )d < B s - [ VsdV^ s , (5.3) 
Jt Jt Jt 



where 



Glials) = ^'[e9 y ( s ^ s + 99 z ( s )r, s + e9 yl ( s )C + e9 zl ( s )r,i}, 



and 



tp(t) = [ T ®'ff v (s)v s + Of v ,(s)v' s ]ds + f T W[fc v (s)v s + 69 v ,(s)v' s ]dB s . 
Jt Jt 

Here we denote, for example, 



E'[0ya)&]= / Ofy(s,u,u')t(s,u>)F{du>') 
Jn 

Jn 



Under (H2), from Theorem 3.1 there exists a unique (£t,r? t ) G S" 2 ([0,T];R) x M 2 (0,T;R) satis- 
fying Q. 
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Lemma 5.2 If we denote by 



e Yf — Yt t E Zf — Z t 

vl = — 6, 4 = — vt- 



Then we have 



r T 

lim sup lE\y £ t \ 2 = 0, lim E / \z £ t \ 2 dt = 0. (5.4) 
£ ^°tefo.Tl Jo 



'te[o,T] 

Proof. First we can express y £ and z £ as 

-dyf = E'[f^t)y! + m)z! + fffltf + f e At)z? + ff(t)]dt 

+V>[g £ (t)y £ t + gt(t)z! + f yl (t)y' t £ + g £ z ,{t)z' £ + g/i(t)]d% 
-zfdWt, 
V e T = 0, 

where, for example, we denote 5 = f,g, Y t ^ = + A(y£ w - Y t>w ), u tjU > = Ut,w + K u t^ ~ 

$y(') = I Sy('i Y -,u),Z. )U ,,u. jU ,,Y. }U i,Z. tU) i,u. tU) i)dX, 
Jo 

and 



TO = [<$(■) - o s y(-m,u + K(-) - e^-u^ + \5 £ yl {-) - 0V(-)]e,u/ 

+[<£(•) - ^'(-)k^ + Wv(-) - &v{-)Ku + WA-) - e s v >(-)]v w . 

Applying Ito's formula to \y £ \ 2 on [t,T], we get 



i-T 



2E^ (y £ ,^[f £ (s)y £ s + f £ z (s)z £ s + + f z ,{s)z'! + /?(«)]> ^ 

/•T 

+E y |E'[^(i)yf + ^( S )^ + ^,( S )yf + ^,( S )4 £ + ff f( S )]| 2 d S . 



From (H2), we have 



E|yf | 2 + E ( T \z £ s \ 2 ds < klE [ T \y £ \ 2 ds + C e , 
Jt Jt 



where k is some suitable constant, C e — > as e — > 0. By Grownwall's inequality, we obtain the 
desired result. □ 
Since u(-) is an optimal control, then 

e -i[J(n £ (-))-J(S(-))]>0. (5.5) 

From this and Lemma 15.21 we have the following 
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Lemma 5.3 Let assumption (H2) hold. Then the following variational inequality holds: 
E / E%(s)£ s +T z (s)ri s +T y >(s)C +t'(s)v' s + Us)v s +T v ,(s)v' s ]ds 



o 



+EE , [/ij / (y ,c,l > ox)eo^ + vftw.VW] > o. 



ifere we denote, for example, 



Proof. From the first result of (|5.4p . we derive 

Jo 

+EE'e- 1 f h y ,(Y 0jU} ,Y 0>U} ,){Y e u , -Y 0>UJ ,)d\ 
Jo 

-> EE , [/ij / (i > 0iW ,y ^)^ + M^^O&w], 

where for example, Yb,u> = ^0,u> + ^(Xq u ~ ^b,w)- Similarly, we have e — > 0, 
£_1 { E ^ E'[/ £ (t) -T(t)]dt| 
-> e/ E'[r a (s)6+T 2 (s)77 s + V(s)^+^K S )^+^(«K+^'( S K]^- 

JO 



Thus dUSD follows. 

Now we consider the adjoint equation: 

pt = ^'h y (%^,%^)+¥rh y ,{%^,Y^) 

+ / F 2 (s,,p s ,q s )ds + / G 2 (s,Ps,g s )dVKs - / q s dJ3 s 
Jo Jo Jo 

where 

F 2 (s,,p s ,q s ) = 1E'[6f y (s)p s + 69 y (s)q s +T y (s)] 

+E*[0> / (s)p* s + 99 y , (s)q* s + T y , (s)], 
G 2 (s, Ps ,q s ) = &[ef z (s)p s + 69 z (s)q s +T z (s)] 

+ E*ief z ,( s )p* s + e9 z ,(s) q * s +T z ,( s )). 

Here we denote by, for example, 

E%/(s)= / y(a,o;*,a;)P(dw*), 

E*[(9> y / (a)?:] = y ef yl (s,u;*,u;) P (s,uj*MdL;*). 
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The adjoint equation (|5.7|) is a special form of (|3.2p with bounded coefficients. Under (H2), it 
follows from Theorem 3.2 that (|5.7j) admits a unique solution (pt,Qt)- 

We define the Hamiltonian function : [0,T] xExRxExRxRxExRxE^Ras 
follows: 



H(t,y 1 ,z 1 ,v 1 ,y 2 ,Z2,v 2 ,p, q) 
= e f (t,uj,uj',yi,zi,vi,y2,Z2,V2)p + g (t,uj,uj',y 1 ,zi,vi,y2,Z2,V2)q 

+l(t,u,u;',y 1 ,zi,v 1 ,y2,Z2,v 2 ). (5.8) 

From variational inequality (|5.6p . we can state the stochastic maximum principle of optimal 
control problem for MF-BDSDEs. 

Theorem 5.1 (Stochastic maximum principle) . Let (Y {■) , Z {■) ,u(-)) be an optimal triple of the 
control problem { (15. ij) . (j5.£j) }. Then Vi> € U, a.e. t € [0, T], a.s. 



\E'H v {t,uj,uj') + ^*H vl {t,uj\uj)} ■ {v-u t ) > 0, 
where for convenience we denote by 
H(t,uj,uj') 

= H(t,u,uj',Yt(oj), Z t (u}),u t (u}),Y t (uj'), Z t (J),u t (J),p t (uj),q t {u)). 

Proof. Applying Ito's formula to {£t,Pt), we obtain 

-E£ Po = E f TS'ilyis^s + Us^s + lyis^ + h'Mlds 
Jo 

-E f F,'[0f v ,(s)v' sPs + 99 v/ (s)v' s q s }ds 
Jo 

-E [ m'[9f v (s)v s p s + 99 v (s)v s q s ]ds. 
Jo 

Then by definition of E* and variational inequality (|5.6p above, we have 

E f E,'[9f vl (s)v' s p s + e9 v r(s)v' s q s +T v >{s)v' s ]ds 







+E / E'[ef v {s)v sPs + 69 v ( s )v s q s + l v (s)v s ]ds > 0. 





From the definition of Hamiltonian function in (|5.8p 

l-T 



E 



o 



E*H V , (t , Y t (oj* ) , Z t (w* ) , % (w* ) , Y t (w) , Z t (w) , % (w) , (w* ) , % (w* ) ) 



(5.9) 



(5.10) 



+E'#„(t, y t ( w ), z t (w), StH, y (u/), Zt(u/),«*K>,M^), 



• v t dt > 0. 
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For Vv G U, F be an arbitrary element of the cr-algebra Ft, set 

u s , s G [0,t), 

, . = I v, s G [*, ,t + s), u eF, 

V ^ S) I u s , s G [i, ,i + e), oj G - F, 

^ u s , s G [t + e,T] , 

we have v(s) G £/ a( t. Since ft satisfies 2* + ft G W a d, then by taking vt = vt — u t , we can rewrite 
above inequality as 

rt+e _ _ 

El F y [M'H v (s,oj,oj') + M*H v ,(s,oj*,u)] ■ (v -u s )ds > 0, 

where is defined in (|5.10p . Differentiating with respect to e at e = gives 

El F [E'fl,;(i,a;,a/) +E*^/(t,o;* J a;)] • (v-u t ) > 0, 
and (|5.9p holds naturally. □ 



6 One mean-field backward LQ problem 

In this section, we are dedicated to apply the previous maximum principle to one backward 
doubly stochastic LQ problem of mean field type. In this case, by supposing h(Yo tU ,, ^o,w') = 

f(s,u),u , Y s ^, Z s ^, ^s,£Jj ) ^s,ui' > Vs,uj' ) 

= Ajy a>w + b^ S)£J + cX* + a 2 iv + s s 2 ^ SiW , + c 2 v S)U/ , 

g(s,LO,u]' 

= Djn, w + ££z s ,w + F.V* + D 2 y S)£ y + E 2 S Z S ^ + F^', 
1(s,uj,uj' , Y S)U }, Z S)U , v SjU1 , Y s ^i , ^ S)a) ' , i> S) o/) 

= ^K 1 ^ + NlZ 2 ^ + fljt,^ + M 2 ^ + N 2 s Zl„, + R 2 s vlJ, 

with, for example, A 1 : [0, T] x f2 2 — > R being bounded, (s,a?,a?') h-> Aj(s, oj, oj') being im- 
measurable (such assumption also hold for the other coefficients), M % ,R l being nonnegative, R l 
being positive, we can write the state equation and the cost functional as 

Jt Jt 
+ y T {[E'F s > s ,^ + B'lDX^' + + FS^d^s,. ~ f Z° a ^w 

(6.1) 
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and 



+^E (£ E,'[M 2 \Y S V >U1 ,\ 2 + A^J 2 + ^I^IV + E'[Q§|r t 



0,w'l J 



For convenience, we write the random coefficients, for example, A 1 (s.oj,oj') as A\ in the above 
and the following part. The Hamiltonian in such setting becomes 

H(s,uj,uj' ,y 1: z 1 ,v 1 ,y 2 , z 2 ,v 2 ,p, q) 
= [A\ yi + B\z x + C\v x + A 2 s y 2 + B 2 s z 2 + C 2 s v 2 ]p 
+ [D] yi + E\z x + FgVi + D 2 y 2 + E 2 z 2 + F 2 v 2 ]q 

+ \[Mly 2 + iV^f + flit,? + M 2 y\ + N 2 z 2 + R 2 s v 2 ]. 

It follows from Theorem 15.11 that 

= E'[C S V + F}q. + Rlu s ] + E* [C 2 ^ + F 2 q* + fl 2 £ s ], 



where, for example, 



E*[fl 2 £ s ] = / A 2 (s,w*,u;)u(,s,u;)P(cL;*), 
E*[C S 2 P :]= / C 2 (s,oj*,oj)p(s,u;*)^(^*), 



and 



+ 



with 



F 2 (s, ,p 8 , q s ) = E'[Alp s + D\q s + M l s ] + ^[A 2 s p* s + D 2 s q* s + M 2 ], 
G 2 (s, p s , q 3 ) = Ef[Bl Pa + E)q a + iV, 1 ] + E* [B 2 p* s + E 2 q* + N 2 }. 



(6.2) 



/ F 2 (s,p s ,q s )ds +/ G 2 (s,p s ,q s )d^ s - / q s dB~ s , (6.3) 
Jo Jo Jo 



Theorem 6.1 Suppose there exists u satisfies h6.S\) , where (p,q) satisfy 116. 3\) , then it must be 
the unique optimal control of above backward LQ problem. 
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Proof. First we have 

J(v) - J(u) 

= \V j\'[Ml{\YlJ - \Y S J 2 ) + N}(\ZlJ - \Z s ^)]ds 
+^E j\'[Rl(\v s , w \ 2 - \u s ^\ 2 ) + M 2 S {\Y^,\ 2 - \Y s ^\ 2 )]ds 
+ 1e j\'[N 2 {\Zl^\ 2 - \Z S ^\ 2 ) + R 2 (\v s ^\ 2 - \u s ^\ 2 )]ds 
+\VM'[QI{\YIJ 2 - \Y^\ 2 ) + Ql(\Y v ^\ 2 ~ \Y y\ 2 )] 

> e / ^[MlY s ^{Y^-Y s ^) + NlZ s ^(Zl^-Z s ^)}ds 
Jo 

+E / E'[fli« 8iW (« a>w - « S)W ) + M 2 y, - y a>w /)]d s 

J o 

+E / E'[iV s 2 Z s ^(^, - Z 8)W ,) + i^« 8)W '(«^ " u a ^)]ds 
Jo 

+EE'[Qjy ^(^ - r 0>w ) + Q^o,u,'(^oV " 
On the other hand, by using Ito formula to p s ,w(Xs,u ~ Y s ,w) on [0,T], we have 

= E / E'[CX" + F, 1 ^]^ - n SjtJ )d S 
Jo 

-E / E'lAfiy^C^ - Y a „) + NlZ^Z^ - Z a , u )]ds 
Jo 

Jo 

-e / e'[m 2 y s ^(y s ^, - y s ^) + jv^v(^ - z S)W 0]<fe. 

JO 
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Thus we have 

J{v) - J{u) 

> E f Ef[Cl PajU + F}q a ^](v ejU -u afi) )d8 
Jo 

+E / Rf[R~u 3 , u (v aiU - u s , u ) + R 2 s u SyLU i(v SyLU > - u s ^)]ds 
Jo 

fT 

+E / E'[C S 2 (> SA / - u s ^)ps^ + F]:(v s ^> - u S)U1 >)q S}U] ]ds 
Jo 

= E/ E'fCgPg^ + F a q St u](v SjU} - u s ,Lu)ds 
Jo 

fT 

+E / [\E'R],]v 8i u(v sjU - u s ,co) + ^*[R%]{v s ^ - u s ,u)]ds 
Jo 

+E / E*[C S 2 P : + F*q*}(v s>w - Us^ds. 
Jo 

Thus by (|6.2p we have J(v) — J(u) > 0, which means u is an optimal control. Since the method 
of proving the result of uniqueness is classical and similar to the case in |13] , we omit it here. 
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